#include "matrix_tool.h"
using namespace arma;
int main(int argc, char const *argv[])
{
	/* code */
	arma::vec A = arma::randu<arma::vec>(5);
	arma::mat X = toeplitz(A,zeros<rowvec>(5)).t();
	X.print("X=");
	arma::mat C = toeplitz_inverse(X);
	C.print("C=");
	X.i().print("X_i=");
	mat D = X*C;
	D.print("D=");
	return 0;
}